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The analysis for the impedance of a vertical electric dipole in the presence of an isotropic 
and homogeneous conducting half-space is presented. Various approaches to the problem 
are then biiefly compared and some numerical results are presented in graphical form. 
The extensions to an anisotropic half -space are also considered. Finally, the dipole is located 
in the space between a homogeneous ground and a sharply bounded ionosphere. It is 
concluded that the presence of the ionosphere has a negligible effect on the impedance of a 
ground-based antenna unless the frequency is less than 1,000 c/s or so. 



1. Introduction 

The primary inotivation for the present investiga- 
tion was the idea that the ionosphere may influence 
the impedance of a ground-based antenna. Nor- 
mally, one would expect the effect to be negligible 
since the separation between the ionosphere and 
ground is large compared with a wavelength for 
medium and higher radio frequencies. However, 
for lower frequencies, particularly those in the 
VLF and ELF regions, the height of the ionosphere 
is comparable with the wavelength. Thus, it would 
appear worthwhile to consider this question. 

A related and more basic problem is the impedance 
of an antenna in the presence of a single boundary 
of separation between two media. Therefore, this 
problem is considered first. 

2. Formulation 

A vertical electric dipole of lengtli ds is located in 
air at height Zq over a homogeneous half-space of 
refractive index N. Tlie magnetic permeabihty of 
the whole space is /x and assumed to be a constant. 
With respect to a cylindrical coordinate system 
(p,</),3), the interface is the plane 2 = 0, the dipole is 
located at z = Zq on the z axis, the half -space ^>0 is 
air with dielectric constant eo, and the half-space 
z<CS) is a homogeneous medium with (complex) 
dielectric constant A^^eo. The situation is illustrated 
in figure la. 

When the current in the dipole varies as J^e'"^^ the 
fields can be found from the solution of a well-known 
boundary value problem [Sommerfeld, 1949]. The 



Hertz vector n has only 
^>0 is given by 



a z component and, for 



n,-ni^)+An3 



(1) 



wh(M"(^ th(^ primaiy influence is 

yciv)_ Jpds exp [-ik\p^+iz—Zoy]^ . 



and the seconchirv inlhuMice is 



(3-.o) 



2\h 



where 



(3) 
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/"= (eoiUo)'<^=27r/wavekM]gth. 



To satisfy the radiation conditions at infinity the 
real parts of u and Ui must be positive for X ranging 
from to infinity.^ 

To carry out a complete calculation of the self- 
impedance of the source dipole requires that the 
electric field in the gap on the connecting transmis- 
sion line be determined from an integral equation 
formulation. In this general case it is necessary to 
consider the finite length of the antenna. However, 
if attention is confined to the change of the imped- 
ance resulting from the presence of the lower half 
space it is permissible to retain the dipole 
approximation. 

B}^ definition the impedance change AZ is related 
to the impedance Z by 



AZ=Z-Zn 



(5) 



1 It may be remarked that if the integration .contour is deformed the real parts 
of u or ui can become negative. This appears to be the source of much confusion 
on the subject. For the present, the contour is the real axis of X and the real parts 
of u are positiv^e. 
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where 



=^]^o= = 



is the free space impedance. It is assumed that Zq 
is either known or can be separately accounted for. 
According to the ''EMF method'' it^follows that 



LZ- 



LEds~\ 



Jo J^ 



(6) 



p-^O 



where thehmits are taken after performing the opera- 
tion 



Ai?.=(F+^,)An.. 
This leads readily to the result 

4xteo'oJo « 



(7) 



(8) 



where q;=22o- 
b}^ writing 



It is convenient to normalize this 
^Z=RJ (9) 



where B.Q = 20¥{dsY is the real part of Zq and is the 
radiation resistance of a dipole in the free space, and 
/ is then a dimensionless parameter which contains 
the essence of the problem. 

3. Asymptotic Approximation 

To obtain an asymptotic approximation to the 
integral for large values of a it is permissible to 
expand R{\) about X=0 in a series of the form 

(10) 

where the prime indicates a differentiation with 
respect to X before setting X=0. As it happens, in 
this particular case, i?'(0), i?'"(0), and all odd- 
ordered differentiated terms vanish. The integrals, 
to contend with, are then of the form 



-f 



\2m+\ 



6-^«(/Xfor m-1,2,3. 



Now 



-r 



— \d\=- 
u a 



(11) 



02) 



and, noting that X^ = «^+i^, it readily follows that 



ika 



This results in the expansion 

\7 / ri^x -ik 



ika 



2! 

/?"'(0) 
4! ^ 



^'+0 



(14) 



and so on. Retaining just the first two terms, this 
may be explicitly written 



/; 

where 






af kajji {kaf' 



i^(0)=^^, 



Fi?"(0)=- 



2 N- 



NN+l 



1 2 

=-rr/^'(0)- 



N' 



(15) 
(16) 

(17) 

(IS) 

(19) 



A question immediately arises as to the validity of 
asymptotic expansions of this kind. In the present 
problem there is a pole where the denominator '>f 
i?(X) is zero. This occurs when 



N\\^~¥)'- + {\^-N'k-'y'=G. 



(20) 



Also, there are branch points at \^=^Nk. To stud}^ 
the influence of these singular points the integral 
P2m is written in the equivalent form 



-I 

Jcq 



X^m^-ua^y^ 



(21) 



where the contour Cq, in the u plane, is sketched in 
figure lb by the solid curve. The poles of the in- 
tegrand, which occur when N^u-\-Ui = 0, are denoted 
by ±Up and the branch points, which occur at u= 
± -x/N^—lk, are denoted ±'^6. The location of the 
singularities in the complex u plane are representa- 
tive for a highly conducting medium where 

iV^!iV|6--/^and !A^^|»1. 

Now the asymptotic evaluation of the integral by a 
steepest descent method amounts to deforming the 
contour Co to contour Ci shown in figure 1. Here, 
(7i is a straight line running parallel to the real axis 
and is defined by the relation u=ik + x where x runs 
from to 00 through real values. 

It is evident from figure lb that no poles or other 
singularities are crossed in the deformation of the 
original contours to the steepest descent contour. 
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Figure la. Dipole over conducting half-space. 



Figure lb. The complex u plane. 



Such is always the case for media in which the real 
part of N^^l. In the case of certain plasma media 
this condition may be violated and tlien the integral 
along Co is asymptotically equal to the integral along 
Ci plus 27ri times the residue of the pole at Up. In 
this study, it is assumed that the real part of iV^>l 
so that the pole contribution does not appear in the 
as3^mptotic evaluation of the integral P2m. 



4. Numerical Integration Approach 

Unfortunately, in many cases of interest, the 
asymptotic series development for the integral is 
poorl.y convergent. Such is the case when |A^| is 
not reasonably large with respect to unit3^ In 
this case, it appears that a numerical integration 
is in order. After some study, the following method 
was employed. A new variable g, defined by \=kg 
is introduced. Then / can be expressed in the form 



T-- 



where 



i' 



Pig) 
io-i) 



I do 



.3 



Pig)=i^R.ikg)g'il+g)-iexp 



(22) 



-ka(g'-~m (23) 



and 



Rikg) 



_NHg'-l)i-(g'-mi 
N'{g'-l)i+{g'-N¥ 



(24) 



To comply with radiation conditions it is necessary 
to choose the square roots so that 



and 



Also 



and 



Lim (g' 
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■^9 


Lim 
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i9'- 
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->(</- 
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(f- 


-N)i- 


^iN 



To simplify the integration further the interval is 
broken into parts such that 



1= 



f 

Jo 



i P(!i)0-u)-'"/g+\ P{g)(ii-^) 



r 



Mg. 
(25) 



These integrals are not particularly suitable for 
numerical integration since there are singularities 
at g=l. However, by further change of variable 
in each integral this objection can be removed. 
In the first integral, we set 



^Vr 



1 



and in the second integral. 

Therefore, 

AZ C"^ r^ 

=/=2 P{l-\-x^)dx-2i P{l-7/)dy (26) 
■^0 Jo Jo 

where the function P is defined by eq (23). These 
integrals are quite suitable for numerical or graphical 
integration. 

5. Compensation Theorem Approach 

Another approach to the determination of A Z is 
based on the compensation theorem from network 
theory. This method has been used extensively in 
the study of antenna ground systems [Wait and Pope, 
1955; Maley and King, 1961]. Using such an 
approach one may write 



where 



AZ. 



AZ-AZ1,_+AZ, (27) 

t=-^2 f H^{pfi)E,{pfi)2Tpdp (28) 

^0 Jo 



where H^ (p,0) is the magnetic field of the dipole on 
a perfectly conducting ground plane at z=0 and 
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Ep(p,0) is the actual tangential electric field on the 
interface z=0. While this formula is exact, it is not 
a solution to the problem since £'p(p,0) is an unknown. 
However, if |A^| is reasonably large compared with 
unity, 

E,(p,0)^-^Hnpfi). (29) 

Inserting this result into eq (28) readily leads to 



AZ 

Rq 



where 



J CO ^-Z 

— 
- ika X 



dx 



(31) 



is the exponential integral. For computational 
purposes it is convenient to use the relation^ 

Ei{-ika) = Ciika) + i[{Tl2)-Si(ka)] (32) 

where Ci and Si are the cosine and sine integrals 
which are extensively tabulated. 

As an interesting check one may use the expansion 



E( 



— IX\_ {—I: 



2! 



3! 



ix) { — ixy {—ixy 



to show, from eq (30), that 

AZ, eg-^^r 1 , i 



Ro 



N 



L(fc 



^,, ^^^^3, terms m^^, etc 



•■] 

(33) 
.]• (34) 



This result may be compared with the first two terms 
of the asymptotic expansion given by^eq (15). 
Noting that 



[f]..r-<'-£)ik"---'"" 



and 



R{0)- 



N-l 



N+l N 



(35) 



(36) 



it immediately follows, from eqs (16) and (17), that 
the asymptotic result is 



N 



lika 



■ay [kaf 
12 



+ 



N 



likaY 



+ 



■■] 
..1. 



-ika 



(37) 



or 



AZt^ — r?— /7 X9 + /7 xo +terms in ,, ,. ; etc. (38) 

which to the order given is identical to the result 
obtained from the method based on the compensa- 
tion theorem. 



6. Discussion of Numerical Results 

Some numerical results based on the above for- 
mulas are now considered. For |A^|=oOj corre- 
sponding to a perfectly conducting plane at 2=0 , 
the impedance increment is given exactly by eq (15). 
The results for this case are shown in figure 2 in the 
form of an Argand plot with the values of ka indi- 
cated on the curves. It is immediately apparent 
that the impedance increment becomes indefinitely 
large as ka becomes small. However, the real 
part of / approaches unity as ka tends to zero.^ 
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Figure 2. The change of impedance of a vertical dipole over a 

horizontal perfectly conducting plane. 

(The results are expressed as a ratio to the free-space radiation resistance.) 



This is physically acceptable and represents the 
doubling of the free-space radiation resistance as 
the dipole approaches the conducting plane. 

Some results for finite values of N are shown in 
figures 3a and 3b for iV=306-^"/* and lOe-'''^^ 
respective!}^. These would correspond to a highly 
conducting half-space where displacement currents 
are negligible. The solid curves are based on using 
two terms in the asymptotic formula for / given by 
eq (15). The dashed curves correspond to just 
using the first term in the same formula (i.e., set 
5=0). The distance between these curves (meas- 
ured on a radius from the origin) gives some idea 
as to the accuracy of the two-term asymptotic 
formula for /. In the case of |A^| = 30, it appears 
that the results should be quite accurate even for 



2 This follows from the relation 



Lim 



i"|«''[-.K'"j)^l=' 
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PiGURE 3a. The change of impedance of a vertical dipole over a 

horizontal imperfectly conducting plane. 

(The results are expressed as a ratio to the free-space radiation resistance.) 

relatively small values of ka. A further check for 
these two cases is afforded by a uiunerical integra- 
tion for ka^?> for each case. This is shown by a 
small circle in both %ure 3a and figure 3b. 

The formula for AZ or / based on the compensa- 
tion theorem is also shown in figures 3a and 3b and 
indicated by broken curves. In hot h cases it falls in 
between the curves for tlie one- and two-term 
asymptotic formulas. Also, it appears to agree 
quite closely with the points obtained from a direct 
numerical integration . 

7. Extension to an Anisotropic Half-Space 

An important extension of the preceding results 
is to a vertical electric dipole over an anisotropic 
half -space. In this case the dielectric constant is in 
the form of a tensor and the problem becomes quite 
complicated. However, if the anisotropic medium 
is a dense ionized medium with a superimposed d-c 
magnetic field parallel to the z axis, symmetry about 
this axis is preserved. In this case, it is permissible 
to employ the Q.L. (quasi-longitudinal) approxi- 
mation. The refractive index p, in this case is 
double-valued and given by 



where 



In tl 



^2^1-i^e^ir (39) 

CO 

tan T=o)jJv, aj,= coo(^^ + a;|,)"2. 
al)ove, coo=angular plasma frequency, 

z^=: collision number, assumed con- 
stant, and 
cojr;,^ angular gyrofrequency. 



Figure 3b. The change of impedance of a vertical dipole over a 
horizontal imperfectly conducting plane. 

(The results are expressed as a ratio to the free-space radiation resistance.) 

It is also assumed that in writing these expressions, 
ijo<^<^v. A good discussion of the validity of the Q.L. 
approximation has been given by Budden [1961]. 

The incremental electric Hertz vector in the air 
space ^^0 is written in tlie following form 



An, 



^ Jods p^A 
47rieocojo 



(X)e7o(Xp)X?/-V/X 



(40) 



which has the same form as cq (8). However, in 
tliis case. 



where 



(pO + Pe)(C'-CoC,) + (pop-l){Co+Ce)C 
{f^O + ^e){C'-Coa)+{pofJ^e + l){Co + Ce)C 

(41) 



c={i-s'y^, s^\/k 

Ce=[i-(s/p,y]^ 

Co=[i-{s/p,y]K 

It can be readily verified that if the d-c magnetic 

A 

field is removed (i.e., t=0) R{X) reduces to Ii(\) 
given by eq (4) since pQ = pe=N. 

The incremental impedance AZ resulting from 
the presence of the anisotropic half-space is tlien 
formally given by 



AZ= 



(dsY 



Jo 



X^?/-V.-"«f/X 



(42) 



4x'ieoco^ 
where again a=2Zo. The two-term asymptotic ex- 
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pansion for this case can be written in the form 

^Z=EJ (43) 

where 

+ 6|P^"(0)|5e^'t^^-s-^"(o)-*^^-^-«l]. (44) 

The complex quantities Ae~''^^ and Be~'^'^f> are de- 
fined b}^ eqs (16) and (17). Using eq (41) for ^(X) 
and carrying out the required differentiations leads 
to 



R(0)= 



MoMf 



-1 



Mo + Me + MoMe + 1 



E'{0)=0 



kUl"{0) = ' 



1 (mo + M6)(2 2 2) 

1 \ Mg Mo/ 

2 Mo + Me + MoMc + l 



(45) 



(46) 



Again, these reduce to the isotropic case when /xo= 

8. Extension to Two-Boundary Case 

Another interesting extension is when the dipole 
is located in a homogeneous space bounded by two 
phine interfaces. The situation is illustrated in fig- 
ure 4. The media beyond the plane interfaces (i.e., 
^<^0 and 2^h) are not specified at the moment. 
However, it is assumed that the fields are everywhere 
symmetrical about the polar axis. 

For the situation described, and for O^^^^o the 
electric Hertz vector has only a z component and 
may be written in the form [Wait, 1960] 



n. 



Jo^.s 



47ri€Q(jo, 



f F(X)Jo(Xp)Xu-V/X (47) 
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Figure 4. Dipole in homogeneous space hounded by two plane 
interfaces. 



where 

mw t^"+-^i(^)^""H^"^^"^^"^+i?2(X)6-^^^--o>] 

^ ^ (?."11+i?i(X)/^(X)6-2"^] ■ ^ ^ 

If 2 and Zq are interchanged hi the above expression, 
the result also holds for the range Zo^Z'^h. The 
coefficients J?i and /?2 may be regarded as reflection 
coefficients for vertically polarized waves incident 
at a (complex) angle sin"^ (X/Z:) at the interfaces z=0 
and z=h, respectively. When the two half-spaces 
are homogeneous media with refractive indices Ni 
and N2, the reflection coefficients have their usual 
Fresnel form given by 



Rj{\)- 



uN)—Uj 



fori=l,2 



(49) 



and u={\^—¥y^ and u^ 



-A:^)^ In anisotropic 
two-bv-two matrices 



cases El and E2 would be 
[Budden, 1961; Wait, I960]. 

A special case worthy of some study is when the 
source dipole is on the bottom interface (i.e., ^o=0) 
and the lower half-space is a perfect conductor (i.e., 
|A7'i|= 00). For this case the incremental impedance 
AZ of the dipole is defined by 

AZ=Z-[Z],= .. 

It follows without difficult V that 



AZ 

2/?o' 



.3 f° 



g^E(kg)e 



-2kh^g 



dg 



[l-E{kg)e-'''^''-'^^g'-^ 



(50) 



where 2i?o is the radiation resistance of the dipole 
located on a perfectly conducting plane. If one 
expands the denominator in a geometric series the 
formal connection between this result and eq (22) 
for the single boundary is evident. This integral 
has not been studied numerically. However, as a 
first approximation one may neglect multiple reflec- 
tions between the two boundaries. This amounts 
to replacing the square bracket term in the denom- 
inator by unity. Then the formula for (AZ/2/?o) is 
the same as that given by the right-hand side of 
eq (15) if a is identified with 2kh. 

9. Concluding Remarks 

The question originally asked has been implicitly 
answered. Namely, the impedance of a ground- 
based antenna is not significantly influenced by the 
presence of the ionosphere unless the reflecting 
heights are somewhat less than a wavelength. For 
example, at a frequency of 15 kc/s it is seen from 
figures 2, 3a, or 3b that the change of impedance 
(relative to the free space radiation resistance) is 
only of the order of one percent. This is hardly 
measurable. However, at 1,000 c/s, ka is of the 
order of 4 and the relative change amounts to about 
15 percent. This is measurable if special precautions 
were made. 
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The important conclusion is that, in transmission 
loss calculations [Norton, 1959], the presence of 
the ionosphere may be neglected insofar as the 
resistance of the antenna is considered. Thus, the 
ionosphere only influences the propagation loss 
[Wait, 1959]. Such an assumption is hnplicitly 
made in most previous studies. 

In a sequel to this paper the effect of raising the 
antenna to a large height will be investigated. 
Actually, the problem has already been formulated 
in this paper. 
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